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Abstract 

We calculate the bound-state energy spectrum of the Dirac Equation 
in a Schwarzschild black hole background using a minimax variational 
method. Our method extends that of Talman £Q to the case of non- 
Hermitian interactions, such as a black hole. The trial function is ex- 
pressed in terms of a basis set that takes into account both the Hermi- 
tian limit of the interaction in the non-relativistic approximation, and 
the general behaviour of the solutions at the origin, the horizon and in- 
finity. Using this trial function an approximation to the full complex 
energy bound-state spectrum is computed. We study the behaviour of 
the method as the coupling constant of the interaction is increased, which 
increases both the relativistic effects and the size of the non-Hermitian 
part of the interaction. Finally we confirm that the method follows the 
expected Hylleraas-Undheim behaviour. 

PACS numbers: 03.65.Pm, 03.65.Db, 03.65. Ge, 04.70.Bw 

1 Introduction 

Variational methods are an important practical technique for determining the 
eigenvalues of Hamiltonian operators. The methods aim to provide an approx- 
imation to the energy eigenvalues by the variation of the Rayleigh quotient 

_ (<W) 

for a continuous set of trial functions <fi — 4>(a). In non-relativistic quantum me- 
chanics the energy spectrum is bounded from bellow and a simple minimization 
of E,), with respect to a leads to an upper bound to the ground state energy. 
Furthermore, if <j> is expressed on a finite-basis expansion, then the eigenvalues 
of H on this basis are approximations to the higher energy eigenvalues. This is 
known as the Hylleraas-Undheim theorem |21 E] . 

For the Dirac equation the picture is more complicated as the spectrum of 
the Dirac Hamiltonian is not bounded from below. That is, the Hamiltonian 
allows states with positive and negative energy. Therefore, a direct minimization 
of the parameter a takes E$ to — oo. This is known as variational collapsc 
13] . Many different techniques have been devised to avoid this problem (for a 
review of the field, see jïïj). Drake and Goldman showed that for a well 
chosen basis the variational approach for the hydrogen atom leads to sensible 
results. They found that the energy eigenvalues split into one positive and 
one negative set. The positive set give the bound-state energies, while the 
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negative energies are smaller than minus the electron's rest mass. Talman pQ 
explained how the finite basis expansion of the large and small parts of the trial 
solution provide the positive and negative eigen-energies respectively. Talman 
also addressed Grant's argument 8 that the variational collapse is caused by 
the incompleteness of the small-component basis. Talman discussed how the 
ground-state energy of a truncated small-component basis is a lower bound for 
the energy of a complete basis. It also turns out that this complete basis energy 
is an upper-bound for the exact energy. Dolbeault et al. jS] have proved that 
equation Q below can be derived from a general variational principle, and have 
stressed that the difhculty in Talman's approach is to determine for which class 
of Hamiltonians it is vàlid. In particular, they have proved that the approach 
is vàlid for Hermitian Hamiltonians with potentials of the form where 
(3 6 (0, 1), which evidently include the Coulomb interaction. 

Recently the bound-state spectrum of an electron in a Schwarzschild black 
hole back ground has been determined ^O]- One of its main features is that 
the energy eigenvalues are complex, since the black hole interaction is non- 
Hermitian. The imaginary part of the energy eigenvalues gives the decay rate 
with time of the corresponding eigenstate. Furthermore, the state spectrum 
resembles that of the hydrogen atom for a small coupling constant. As the 
coupling constant is reduced, the imaginary part of the energy tends to zero as 
the interaction tends to one of Coulomb type. This offers a starting point to 
cxplore the variational methods for non-Hcrmitian Hamiltonians. In this parer 
we thus implement a minimax mcthod based on Talman's approach for the black 
hole interaction. This method allows us to calculate the full complex energy 
eigenvalues and explore its behaviour with an increasing coupling constant. The 
approach developed here should have applications in a range of fields where 
one is forced to deal with non-Hermitian interactions, including the important 
subject of open systems. 



2 Talman's minimax method 

Talman's variational formulation for the energy eigenvalues follows from the de- 
composition of the Dirac wavefunction into its large and small parts. Following 
Grant jS], the two-component radial Dirac equation can be written as 

VJ - U*P + è-A -1 + ^H l{ r)) [f) ~ me* [f) (2) 

where H±(r) is the interaction term. The ground state energy can be obtained 
from the variation of the correct eigenfunction 0, determined by / and <?; 



E 



maxj 



,H0) 



where <fi is the two-component radial spinor 

9{r) 
f(r) 



(3) 



(4) 



If g is an arbitrary function, the maximum over / yields an upper-bound to the 
exact ground energy. The trial function <f> can be expanded in a finite function 
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basis in the form 

After a substitution of this function into equation the minimax approxima- 
tion is reduccd to minimize E for each V and maximize it for each a J . If wc 
then consider vectors of the form 

a = ai, ...On)*, (6) 

then the minimax variation of E leads to a generalized eigenvalue problcm. This 
can be written 

Here h is a (2n x 2ri) matrix 

h = 

where the (n x n) matrices are 

hlly 
hl2y 
h21y 

h22y 

The matrix O is also a (2n x 2n) matrix given by the overlap integrals of the gi 
and fj. Its hrst (n x n) sub- matrix is 

/•OO 

Olly = / gfgjSdr, (13) 
Jo 

with 022 containing the equivalent fj integrals, and all other terms vanishing. 

Talman shows that for the Coulomb interaction the positive eigenvalues in 
equation J7J tend to the first energy eigenvalues as n increases, following a 
Hylleraas-Undheim behaviour. 



ha = EOa. (7) 

(8) 



(hllnxn hl2 nX ra 
h21 IlX n h22 nX n 



given by 



/>oo 

/ 9i Hu gj r 2 dr 
Jo 

poc 

/ gtH 12 f jr 2 dr 
Jo 

/ /* H 2 i gj r 2 dr 
Jo 



(9) 
(10) 

(11) 



f?H 22 f jr 'dr. (12) 



3 Black hole interaction 

The interaction Hi for a black hole is given by ^I] 



1 TT . . . h h 2GMm 3 , 

— Hx{r )=%— J — (d r + — ), 14 

mc z mc V mer en 4r 

where the dimensional constants have been arranged into a convenient form. If 
we measure the radial distance in units of the Compton wavelength (mc/h) we 
are left with the dimensionless coupling constant 

GMm . . 
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which gives the magnitude of thc intcraction. Using the interaction (|14|) . we can 
now write the system of differential equations (0) in a more convenient form. 
With the change of variables 



/ 

w\ = — 

r 

the equations can be written as 



where C is the matrix 



d r vo\ 

d r ZÜ2 



and 



= C(r) 
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= — , 

r 



( 



(16) 



(17) 



C(r) = 



2a 



E) 



i(-k 




E) 

'2ar{\ - E) 

(18) 

and we measure r and E in units of (mc/h) and (mc 2 ) respectively. To simplify 
some of the later expressions we set the Compton wavelength mc/h equal to 1. 

Special attention must be given to the singularities of equation lfT?|) . The 
equation has three singular points. Two of them are regular singularities sit- 
uated at the origin and at the horizon, r = 2a. The third one is at infinity 
and is an irregular singularity of rank 2. Expansions of the solutions around 
the singular points provide the boundary conditions for the differential equa- 
tions. The physical solutions are then selected from the correct boundaries at 
the singularities. Gathering the complete behaviour of the bound-state solution 
around the singular points, we obtain the general form 



ip(r) = R(r)e 



-Ví-E 2 r+2iV2arE , 



-3/4 



(19) 



where R{r) is an analytic function everywhere. The exponential factor and 
the first term in r take care of the singularity at infinity |12|. The r -3 / 4 term 
corresponds to the boundary at the origin. The presence of this term is crucial 
as it accounts for the non-Hcrmiticity of the Hamiltonian JI] . This can be seen 
more explicitly from the calculation of the matrix elements of hll and h22. 
We find, for example, that 



hll. 



gf g 3 r 2 dr + V2a~ gf r^Hd r (r 3 / 4 9j )] dr 



g* gir 2 dr 



+ iV2a~r 3 ' 2 gïg 3 \l 



g*r 3 /Hd r (r^gt)dr 



hll*,-/ + i 



9t 9 j\o ■ 



(20) 



The second and third terms of the second equation follow from an integration 
by parts. A similar calculation holds for h22y. As a consequence, if there is at 
least one pair gi and g 3 where both tend to the origin as r~ 3 ' 4 , thc Hamiltonian 
is non-hermitian and the energy complex. The imaginary part of the energy is 
always negative and represents the decay rate of the wavefunction with time. 
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4 Minimax for the black hole interaction 



Despite the fact that the interaction Hamiltonian is non-Hermitian, we can 
still apply the minimax approach, if we select a suitable function basis. The 
choice is based on the assumption that we achieve the non-relativistic limit for 
low a, in analogy with the Hydrogen atom. Lasenby et al. 10 show that 
the non-relativistic limits of the black hole and Coulomb interactions coincide. 
Therefore, we expect that the energy spectrum for the black hole resembles the 
Hydrogen spectrum for a low coupling constant. Restricting ourselves to the 
low coupling regime, the task is then to choose a finite function set that provides 
the Hydrogen properties and, at the same time, has the correct boundaries at 
the singularities. 

From equation (|19fl we see that the bound state solutions of the Dirac equa- 
tion must go as r -3 / 4 for r — > 0, and as e - br + l PVr f or r _> OQm Taking these 
into account we define a set of discrete basis functions by 

g = e lVïï ^e~ br r-i gj = é v^F e -&V-i ( 2 i) 

and 

/„ = e iV2^ e -«r r -| j. . = e iV25F e -«yí-i (22) 

where j = l...n, and a and b are additional (nonlinear) parameters. We limit 
ourselves for the case where a and b are real. This corresponds to the condition 
^s(E) << $l(E), expected for the non-relativistic limit. Note that we have 
increased the order of the approximation to n + 1 by including the extra terms 
go and /q. From equation (|20[1 . the matrix h will non-hermitian if we introduce 
go and fo, and we find the the non-hermitian terms of h from the quantitics 

^9Wr)g r^r = _ ab _ ^ 
Jo 9o 9o r 2 dr \ ir 

and 



lo 00 fo for 2 dr V t 

For the remaining terms down the diagonal of h we find that 

í°° 

g*H 1 (r)g 3 r 2 dr = - J g*(-jg jr 2 dr (25) 

with an identical result holding for fj . We can see here the Coulomb interaction 
explicitly appearing in the diagonal of h. 

Leaving out the go and fo terms in the expansion, the minimax method 
reduces to a real energy approximation. Thus, applying the minimax method 
to a pair of functions gj and fj with (j ^ 0) allows us to fix the vàlues of a 
and b. With these vàlues, the full basis can be determined, including the case 
j = 0. Our central assumption now is that equation (Q is remains vàlid, even 
when using the full basis set with h non-Hermitian. The generalized eigenvalue 
problem of equation (JJJ allows us then to calculatc complex eigen-energies. The 
validity of our assumption rests in part on the fact that the results obtaincd 
are in good agreement with the results of othcr, more accurate techniques. 
Since equation {JJ can be derived from a variational principle of a Hermitian 
Hamiltonian, its validity for a non-Hermitian Hamiltonian suggests the existence 
of a variational principle for such a case. 
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Figure 1: Rayleigh quotient as function of the non-linear parameters a,b. a = 
0.1, k = -1; minimax at: a = 8.9186905 x 10" 2 , b = 7.8238289 x 10~ 2 . 

5 Results 

We start with thc normalizcd basis function 



which sets j = 1 in equations l|21|) and l|22|) . As discussed, the Rayleigh quotient 
for this function will be real. So the solution of equation Q is given by the 
saddle point of as a function of a and ò. An example of this saddle point is 
shown in figure ^ From equation l)25|l we see that the same result would have 
been obtained for this function basis had we used the the Coulomb interaction 
instead. 

The vàlues of a and b determine a unique function set specified by equations 
(I21|l and 122|) . A suitable trial function expanded on this function set is then 



With this trial function, the eigenvalues of equation JJJ) can be obtained. Since 
the function with power —3/4 in equation Ij27(l gives rise to the non-Hermitian 
terms in h, the eigenvalues are complex. The approximation to the ground 
statc energy is then obtained from the eigenvalue with minimum positive real 
part. As the order of the expansion is increased, we expect the approximation 
to improvc its accuracy and, in particular, tend to a limit. That this is indeed 
the case is confirmed in figure [3 

The minimax procedure offers only an approximation to the exact bound- 
state energies. The accuracy of the method depends on the chosen trial function. 
We consider two extra trial functions and compare their results with a numerical 
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Figure 2: Convergence to a single value of the approximation to the real part 
of the energy ^St(E). For vàlues a = 0.1, k = — 1, the energy E/(mc 2 ) — 
0.9947208 - i 2.0498160 x 10~ 5 is obtained. 

shooting method that provides accurate vàlues of the ground-state energy 
The shooting method finds the complex energies vàlues for which the integrated 
solutions from the boundaries at the horizon and at infinity coincide. For the 
first trial function we consider the basis functions increasing on 1/4- rs powers 
of r. Thcn the trial function can bc written as 



equations l|27(l and (|28|l . The new trial function increases in 1/4 powers for the 
negative powers of r and in integer vàlues for positive powers. This is 



For a = 0.1, k = — 1, we find the ground-state energies given in table^ (Note 
that the present conventions have k with opposite sign to the convention for n 
in ^ü], so k = — 1 is the ground state.) The table shows that the approximation 
to the real part of the complex energy given by the second trial function (n = 16) 
is extremely good, with the imaginary part also obtained quite accurately. The 
third trial function (n = 18) also gets the real part quite accurately, but the 
decay rate approximation is poor. 

The trial function 02 is then more suitable for the study of the behaviour 
of the method with an increasing coupling constant. This is shown in figures [3] 




(28) 



For the second trial function we consider a combination of the trial functions in 




(29) 
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E/(mc 2 ) 

Trial function - minimax <\> x 0.9947208 - i 2.0498160 x 10" 5 
Trial function - minimax (f> 2 0.9946883 - i 2.7855588 x 10~ 5 
Trial function - minimax (f> 3 0.9946858 - i 3.0659648 x 10~ 5 
Shooting method 0.9946882 - i 2.7870824 x 10~ 5 

Table 1: Ground-state energy. The energies are obtained from the minimax, 
using diffcrent trial functions, and from the shooting method. 



and 01 where we can see that the minimax method gives good results for low 
couplings. Around the point of the last stable circular orbit a — 0.28 [TU] 
the minimax loses accuracy. This gives numerical indication that the region 
where the minimax works at its best is where classical circular orbits are still 
admissible. 

The minimax method also provides approximations for the excited states. 
These are obtained from the other real positive eigenvalues. Using the trial 
function 03 a neat Hylleraas-Undheim behaviour is obtained for 3t(E). This 
is shown in figures and [ÏÏJ where the real eigenstates tend to different energy 
levels as the order increases. In contrast to the hydrogen case, we have not only 
positive energy but also negative energy states. A selection of the vàlues for the 
bound energies obtained in these figures are given in table [21 where the first two 
eigenvalues of positive and negative energies are compared for the cases k = — 1 
and k = 1, with a = 0.1. 

From this table we can extract a number of features. There is a spurious 
root for the first negative state with k = 1. All the other energies can be 
checked in a numerical integration routine and confirmed to give bound states. 
The solution for this spurious root does not show any structure other than 
an exponential divergence from the origin. The presence of spurious roots in 
variational methods for the Dirac equation is not new. In the past these roots 
have appeared from the study of the Coulomb interaction. In previous works, 
Drake and Goldman |üj found the presence of a spurious root for the first 
eigenvalue of the case k = 1. In his original work 1 Talman showed how, for 
a particular case, the minimax method removes the spurious root found in [Tj- 
Note that this kind of root is not present in table so the minimax procedure 
has certainly excluded this case and given a sensible value for the first positive 
eigenvalue with k = 1. The spurious root is found, however, for the first negative 
eigenvalue of k = 1. This lies outside the scope of Talman's arguments, since 
negative energy states of this type are not present in Coulomb interactions. 

Another feature of table is given by the comparison of the energies of 





Positive Energy/mc 2 


Negative Energy/mc 2 


fc = -l l st 

k = l l st 

2^nd 


0.994686 - i 3.06596 x 10" 5 
0.998698 -i 3.83240 x 10~ 6 
0.998731 -i 2.00712 x 10~ 7 
0.999438 -i 5.22900 x 10~ 8 


-0.998730 - i 1.53600 x 10" 8 
-0.999437- i 4.74152 x 10~ 7 
-0.985823- i 1.55033 x 10~ 2 
-0.994978 - i 1.75747 x 10~ 3 



Table 2: Bound- state Energies for excited states 
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Figure 3: Real part of the energy as a function of the black hole mass. The real 
part of the energy given by the minimax (light huc) and by the shooting energy 
(dark hue) are compared for the ground state with k = — 1. 
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Figure 4: Decay rate as a function of the black hole mass. The decay rates 
given by the minimax (light hue) and the shooting (dark hue) corresponding to 
figure 13 
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Figure 5: Hylleraas-Undheim behaviour of the positive real part of the eigen- 
energies, given the vàlues a = 0.1, A; = — 1 and the trial function (|29(l . 
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Figure 6: Hylleraas-Undheim of the negative real part of the eig en- energies, 
given the vàlues a = 0.1, k = — 1 and the trial function (|29ü . 
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the states 2S x / 2 (2 nd ,k = -1) and 2P 1/2 (l st ,k = 1). For the low value of a 
used here, we see that the states have comparable energy, as is the case for the 
Coulomb interaction. At higher vàlues of a for the gravitational interaction this 
degeneracy of the energy in the sign of k is lifted [TÜ] , 

A further property we can extract from table El is the presence of charge 
conjugation symmetry For a given state we should find another with opposite 
sign for both the real part of the energy and k, but with the same decay rate. 
This can be seen from the comparison between the real part of the state (l sí , k = 
— 1) of positive energy and the state (2 nd , k = 1) of negative energy, as well as 
the comparison between the states of positive energy (l sí , k = 1), (2 nd , k = 1) 
and the states of negative energy (l 3t ,k = —1), (2 nd ,k = —1) respectively. It 
is a bit disappointing that the imaginary energies of corresponding states have 
not the same value, as is required by the charge conjugation symmetry. This 
suggest that the minimax method will need to be further refined to find a more 
precise decay rates for the excited states. 

6 Conclusions 

The minimax method has proven to be a useful tool for analyzing the electron 
bound state energy spectrum in a black hole background. We have seen how 
the method works well not only for the small coupling limit but also sheds some 
light on what happens beyond the last classically stable circular orbit. The 
validity of the method suggests that the complex energy spectrum of the black 
hole can be obtained from a general variational principle. 

The minimax method can only give an approximation to the exact energy. 
As we have seen this depends critically on the trial function used. Furthermore, 
we found that the method ceases to give reliable results around a ~ 0.4, which 
may limit the applicability of the minimax method. One reason for this failure 
may be that we have assumed that the imaginary part of the energy is small as 
compared to the real part. This is no longer the case around a ~ 0.4 for the 15* 
state. However, we may expect that the minimax method will continue to give 
accurate results for higher angular momentum states beyond a ~ 0.4. 

A further improvement of the minimax method can be given by the relax- 
ation of the condition on the imaginary energy. As a consequence, the nonlinear 
parameters a and b could be considered as complex numbers. It then would be 
interesting to see whether this gives any improvements in either the determina- 
tion of the energy for higher vàlues of the black hole mass, or to the accuracy of 
the decay rates for excited states. Extending this method to larger vàlues of a is 
a particularly significant problcm, as the shooting method runs into numerical 
difiiculties beyond a of around 6 , whereas vàlues of astrophysical interest start 
at around 10 15 . 
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